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This  briefing  contains  information  up  to: 


Isomers  of  bis(pentazolyl)iron(ll) 
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HEDM  Program  Objective 


Breaking  the 
performance  barrier 


Identifying  and  developing  advanced  chemical 
propellants  for  rocket  propulsion  applications 


•Isp  is  the  major  metric  of  a  propellant’s  performance 
•Density  can  also  be  a  significant  contributor 
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Why  We  Are  Doing  It 
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10%  increase 
in  both 


05 


Propellants  Program 
General  Ap p roach 


Employ  a  synergic  blend  of  experimental  (synthesis  and  physical)  and 
computational  techniques  derived  from  the  disciplines  of  chemistry  and  physics 
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Polynitrogen  Project 


Discover ,  synthesize ,  characterize,  and  scale-up 
novel,  highly  energetic  polynitrogen  allotropes 

Modeling  and_  simulation  guides  the  experimental  program: 

♦  Determines  which  molecules  should  exist  and  how  energetic  they  are 

♦  Gives  information  on  howto  synthesize  promising  molecules 

♦  Provides  critical  data  for  identification  and  characterization  of  new  molecules 
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Performance  of 
Polynitrogen  Monopropellants 


The  performance  of  polynitrogens  as  monopropellants  would  dwarf 
that  of  hydrazine,  and  would  greatly  exceed  even  bipropellants 
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The  Search  for  New  Polynitrogens 


•  All  polynitrogens  are  unstable  with  respect  to  N2  molecules 


•  Their  activation  energy  for  N2  elimination  is  largely  determined  by  the 
weakest  bond  in  the  compound 


•  Their  metastability  is  enhanced  if  suitable  resonance  structures  exist: 

[l  N-N=N|]  4 - >  [1g=N=N^]_  < - >  [|N=N-N|]~ 

•  The  double-bond  character  of  the  N — N  bonds  in  the  azide  anion 
explains  its  exceptional  stability 


•  How  can  this  stabilization  effect  be  used  to  our  advantage  in  preparing 
new  compounds? 


This  page  is:  DISTRIBUTION  A.  Approved  for  public  release;  distribution  unlimited. 


Pentazolate  (N5')? 
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*  Substituted  pentazoles  R-N5  have  been  known  for  decades  (R=aryl) 

*  Cyclic  N5~  is  aromatic 

*  Conversion  of  the  diazonium  salt,  RN2+,  to  the  substituted  pentazole  ring  R- 
N5  by  the  reaction  with  azide  ion,  N3~,  has  been  demonstrated  many  years 
ago  by  Ugi  and  Huisgen. 

*  N5-  has  been  recently  detected  in  the  gas  phase  for  the  first  time,  using 
collisional  fragmentation  (electrospray  ion  mass  spectroscopy). 

*  Can  a  chemical  route  to  N5~  be  found? 
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New  Polynitrogen  Anions 
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Pentazole  anion  (N5*) 

•  Theoretical  calculations  show  that  this 
anion  has  a  28  kcal/mole  activation  energy 
barrier  for  decomposition  and  its 
decomposition  to  N3  and  N2  is  only  1 1 
kcal/mol  exothermic 

•  Aryl  substituted  pentazoles  can  be  isolated 
as  stable  compounds  only  if  stored  at  low 
temperatures.  In  methanol ,  these 
compounds  rapidly  decompose  at  room 
temperature  to  form  aryl  azides  and  N2  gas 
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◄LL# 


:k  Transition  Metal  Complexes  of 


[NJ? 


[N5]*  is  isoelectronic  with  cylcopentadienyl  anion  [C5H5] 
(cp) 


cp  ligands  readily  bind  to  many  transition  &  main  group 
metals  to  form  “sandwich”  complexes 

Ferrocene  (cp-Fe-cp)  -  first  observed  in  1951  (T.J.  Kealy,  P.L.  Pauson 
Nature,  168,  1039(1951);  S.A.  Miller,  J.A.  Tebboth,  J.F.  Tremaine  J.  Chem.  Soc. 
632  (1952).) 

Chemistry  of  metallocenes  has  been  thoroughly  studied 


Insights  into  possible  reaction  pathways  for  isoelectronic  N5-TM-N5 
compounds? 
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[N5]'  is  isoelectronic  with  aromatic  cylcopentadienyl  anion  [C5H5]' 


cp  ligands  can  coordinate  to  metal  center  via  tt  electrons  only  (r|5) 

[N5]‘  can,  in  principle,  coordinate  via  tt  electrons  or  a  lone  pairs 

tt  coordination:  n5 
a  coordination:  n1  or  r\2 
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^  Quantum  Chemical  calculations  of 


N5-Fe-N5 


Determine  possible  bonding  modes  and  their  relative 
energies 


Determine  possible  spin  states  and  their  relative  energies 

Singlet,  triplet,  quintet  electronic  states. 
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4  1  Quantum  Chemical  calculations  of 

Ng-Fe-Ng 

Results  of  initial,  modest-level  calculations  (B3LYP/SBK+(d)  and 

MP2/SBK+(d)) 

Singlet  electronic  states,  relative  energies  in  kcal/mol 


n5-n5 

B3LYP 

21-5  (D5d) 

MP2 

3-6  (D5h) 

n2-n2 

10.1  (D2d) 

4.6  (D2h) 

0.0  (C2) 

10.0(D2d) 

0-5(D2h) 

0.0  (C2) 

rV-n1 

9.2  (C2v) 

n/a 

rV-n2 

n/a 

5.8  (Cs) 

n1-'!5 

n/a 

n/a 

n2-n5 

13.2  (Cs) 

9.5  (Cs) 
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Quantum  Chemical  calculations  of 

N5-Fe-N5 


B3LYP  MP2 

10.1  (D2d)  10.0  (D2d) 


B3LYP  MP2 


13.2  (Cs)  9.5  (Cs) 


B3LYP 

0.0  (C2) 


MP2 

0.0  (C2) 


MP2 

3-6  (D5h) 


B3LYP  MP2 

4.6  (D2h)  0.5(D2h) 


B3LYP 

21-5  (D5d) 


5-8  (Cs) 
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◄LL# 


^  Quantum  Chemical  calculations  of 


N5-Fe-N5 


B3LYP,  triplet  electronic  states,  relative  energies  in  kcal/mol 


n2-n2 

n2-n5 
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■12.1  (Cs) 
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Total  (au)  and  Relative  (kcal/mol)  Energies  of 
N5-Fe-N5  stationary  points 


B3LYP(5) 

MP2 

CCSD(T)//MP2 

T1 

n5-n5 

D5h 

-221 .897766[1  ](21 .5) 

-220.722672[0](  3.6) 

-220.817543(12.2) 

0.0435 

D5d 

.897778[0](21.5) 

.721 345[1](  4.4) 

-220.818049(11.9) 

0.0432 

02^2 

D2h 

-221 .91 4555[1  ](1 0.9) 

-220.71 1 445[1  ](10.6) 

-220.83341 8(  2.2) 

0.0253 

C2h 

.924379[1](  4.8) 

.721 884[1  ]( 4.1) 

n/a 

n/a 

D2h* 

.924679[0](  4.6) 

.727578[0](  0.5) 

-220.831 820(  3.2) 

0.0218 

D2d 

.91 5901  [0](1 0.1 ) 

.712431  [0](1 0.0) 

.833717(2.0) 

0.0254 

C2 

.931 975[0](  0.0) 

.728409[0](  0.0) 

.836247(  0.4) 

0.0296 

D2d(triplet) 

.91 6963[0](  9.4) 

D2h(triplet) 

.97871 5[0](-29.3) 

D2h*(triplet) 

.981686[0](-31.2) 

nl-ni 

D2h 

-221.869702[2](39.1) 

-220.676205[2](32.8) 

-220.788778(30.2) 

0.0277 

D2h* 

.906936[1  ](1 5.7) 

.7011 73[1  ](1 7.1 ) 

-220.812963(15.1) 

0.0257 

D2h** 

.896858[2](22.0) 

.707632[1  ](1 3.0) 

-220.809089(17.5) 

0.0188 

D2d 

.847741  [2](52.9) 

.643839[2](53.1) 

-220.789710(29.7) 

0.0228 

C2v 

.91 7363[0](  9.2) 

.711 637[1  ](1 0.5) 

-220.817459(12.2) 

0.0246 

Cs 

n/a 

.692082[1](22.8) 

n/a 

n/a 

D2h(triplet) 

.952592[0](-12.9) 

c2V(triplet) 

.945762[0](-8.7) 
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Total  (au)  and  Relative  (kcal/mol)  Energies  of 
N5-Fe-N5  stationary  points 


B3LYP(5) 

MP2 

CCSD(T)//MP2 

T1 

nl-n2 

C2v 

-221 .898673[2](20.9) 

-220.7021 1 3[2](1 6.5) 

-220.807136(18.7) 

0.0223 

C2v' 

.914327[1](11.1) 

.713738[1](  9.2) 

-220.820661(10.2) 

0.0230 

Cs 

.91 9090[1  ](  8.1) 

.7191 27[0](  5.8) 

-220.823692(  8.3) 

0.0236 

Cs(triplet) 

.951 282[0](-1 2.1 ) 

ni-n5 

Cs 

n/a 

-220.688408(1  ](25.1 ) 

-220.801039(22.5) 

0.0402 

n2-n5 

Cs 

Cs‘ 

-221.910857[1](13.3) 

.91 0869[0](1 3.2) 

-220.71 3209[0](  9.5) 

-220.836962(  0.0) 

0.0451 

Cl  (triplet) 

.943231  [0](-7.1) 
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Occupation  Restricted  Multiple 
Active  Space  (ORMAS) 


ORMAS:  flexible  CAS  method  for  reducing  size  of  full  Cl  step, 
similar  to  generalized  CASSCF. 

Active  space  is  product  of  three  sub-spaces: 

(6e,5o)  [N5]-  tt+  tt*  (left) 

(6e,6o)  [Fe]2+  3d+4s  (middle) 

(6e,5o)  [N5]-  tt+  tt*  (right) 

SBK+(d)  effective  core  potential  and  valence-only  basis  set 
used  throughout. 
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ORMAS  Results 
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-219.351299 
(0.0  kcal/mol) 


Natural  Orbital  Occupation  Numbers 


N5  (left) 

Fe 

N5  (right) 

tt  1.983 

3d  1.993 

tt  1.976 

TT  1.938 

3d  1.000 

tt  1.935 

TT  1.916 

3d  1.000 

TT  1.930 

TT*  0.086 

3d  1.000 

TT*  0.083 

TT*  0.077 

3d  1.000 

TT*  0.077 

4s  0.007 
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ORMAS  Results 
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C2v 

n2-n2  r^) 


-219.182241 
(+106  kcal/mol) 


Natural  Orbital  Occupation  Numbers 


N5  (left) 

Fe 

N5  (right) 

tt  1.976 

3d  1.989 

TT  1.996 

tt  1.936 

3d  1.000 

tt  1.855 

TT  1.931 

3d  1.000 

tt  1.027 

TT*  0.082 

3d  1.000 

TT*  0.974 

TT*  0.076 

3d  1.000 

TT*  0.148 

4s  0.011 

Distribution  A:  Approved  for  public  release;  distriubtion  unlimited. 


ORMAS  Results 
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^2h 

n2-n2  (3B2g) 


-219.260755 
(+57  kcal/mol) 


Natural  Orbital  Occupation  Numbers 


N5  (left) 

Fe 

N5  (right) 

tt  1.977 

3d  1.963 

tt  1.977 

TT  1.935 

3d  1.931 

TT  1.935 

TT  1.930 

3d  1.035 

TT  1.930 

TT*  0.082 

3d  0.982 

TT*  0.082 

TT*  0.077 

3d  0.084 

TT*  0.077 

4s  0.006 
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Natural  Orbital  Occupation  Numbers 


N5  (left) 

Fe 

N5  (right) 

tt  1.972 

3d  1.972 

tt  1.972 

TT  1.951 

3d  1.900 

TT  1.951 

tt  1.949 

3d  1.024 

tt  1.949 

TT*  0.066 

3d  0.985 

TT*  0.066 

TT*  0.064 

3d  0.114 

TT*  0.064 

4s  0.005 
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ORMAS  Results 
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Cs 

n5-n5  ('A1) 


2.221(4) 


-219.127764(7) 
(+140  kcal/mol) 


Natural  Orbital  Occupation  Numbers 


N5  (left) 

Fe 

N5  (right) 

tt  1.972 

3d  1.996 

tt  1.972 

tt  1.945 

3d  1.834 

tt  1.945 

tt  1.945 

3d  1.824 

tt  1.945 

tt*  0.070 

3d  0.171 

TT*  0.070 

TT*  0.068 

3d  0.171 

TT*  0.068 

4s  0.004 
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Conclusions 
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•  At  the  B3LYP  and  MP2  levels  of  theory,  there  are  six  distinct  stable  singlet 
structures.  The  most  stable  structure  is  a  twisted  r|2-r|2  isomer  with  C2 
symmetry. 


•At  the  B3LYP  level,  there  are  seven  distinct  stable  triplet  structures.  The  most 
stable  structure  is  a  planar  n2-n2  isomer  with  D2h  symmetry. 


•At  the  B3LYP  level,  all  but  one  of  triplet  minima  are  more  stable  than  the 
lowest  energy  singlet  structure.  The  most  stable  triplet  is  a  planar  n2-n2  isomer 
with  D2h  symmetry. 


•ORMAS  results  indicate  significant  multiconfigurational  character  in  N5-Fe-N5 
complexes  -  B3LYP  and  MP2  results  may  not  be  reliable. 
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